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We show that the phase appearing in neutrino flavor oscillation formulae has a geometric and 
topological contribution. We identify a topological phase appearing in the two flavor neutrino 
oscillation formula using Pancharatnam's prescription of quantum collapses between nonorthogonal 
states. Such quantum collapses appear naturally in the expression for appearance and survival 
probabilities of neutrinos. Our analysis applies to neutrinos propagating in vacuum or through 
matter. For the minimal case of two flavors with CP conservation, our study shows for the first 
time that there is a geometric interpretation of the neutrino oscillation formulae for the detection 
probability of neutrino species. 
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I. INTRODUCTION 

The phenomenon of neutrino flavor oscillation results 
from the phase difference acquired by the mass eigen- 
states due to their time evolution while propagating in 
vacuum or in matter. The observation of neutrino fla- 
vor oscillations in solar, atmospheric, reactor, and accel- 
erator experiments reveal the remarkable fact that the 
neutrinos exhibit sustained quantum coherence even over 
astrophysical length scales P, 0]. It is then natural to 
ask what we can learn about neutrinos from these co- 
herent phases. Here, we address the issue of geometric 
and topological phases involved in the physics of neutrino 
oscillations. 

On the theoretical front, it is well known that the 
phenomenon of neutrino oscillations cannot be accom- 
modated within the standard model (SM) of particle 
physics. Therefore, the experimental observation of neu- 
trino oscillations provides a concrete evidence for the re- 
quirement of physics beyond the SM and neutrinos have 
been an intensive area of research in the past several 
years. 

The study of geometric phases in the context of neu- 
trino oscillations has been carried out in t he p ast by sev- 
eral authors @, i, n, i, 0, s, s m m m, m m m m, 

UtI [Tc| . but none of the papers seem to provide a uni- 
fied perspective on the problem taking into account the 
different avatars of geometric phase. It is worthwhile to 
stress here that one needs to be cautious while interpret- 
ing claims in the literature as they crucially depend on 
which version of the geometric phase one is dealing with. 
We will first summarize the related literature and then 
focus on the specific question that we address in this 
paper. We mostly restrict our attention to the case of 
two neutrino flavors and the CP (CP stands for charge- 
conjugation and parity) conserving situation, which is the 
minimal scenario for studying the physics of oscillations. 
We find, in contrast to earlier studies of this problem 
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that the geometric phase appears even in this minimal 
context. 

Let us first review the papers that are connected to 
Berry's cyclic adiabatic phase. Berry studied phases 
that appear when the Hamiltonian of a quantum sys- 
tem depends on parameters that are varied slowly and 
cyclically. Nakagawa @ followed this work by an ele- 
gant paper in which he pointed out that the geometric 
phase could also arise in systems where adiabatic theo- 
rem did not hold. The key point made by Nakagawa was 
that while for existence of geometric phases, adiabatic 
condition was not necessary (this was also independently 
pointed out by Aharonov and Anandan 20]), the adia- 
batic theorem itself could be most easily understood in 
terms of geometric arguments. As an application of his 
general formalism, Nakagawa considered two flavor neu- 
trino oscillations in matter. He concluded that the Berry 
phase played no role in this situation. The topological 
phase in the two flavor neutrino case, which is the cen- 
tral result of the present paper was missed in his work 
because he restricted himself to a limited region in the pa- 
rameter (ray) space and did not consider generalizations 
of the geometric phase that allow for quantum collapse. 

Subsequent work on Berry's geometric phase and neu- 
trinos exploited the spin degree of freedom of neutrinos 
and its interaction with the transverse magnetic field 
leading to geometric effects and spin flip. Since at that 
time, spin precession was a plausible solution to the solar- 
neutrino problem, there is a body of work by several au- 
thors on the subject of geometric phase effects in this 
context, both in the absence and presence of matter and 
mass-splitting terms [4|, la, |g, |7|, [8j, l9j, llQ|> 111! • However , in 
the present scenario, spin flavor precession is disfavored 
as the leading solution to the solar-neutrino problem at 
99.86 % C.L. [2l[, which makes it phenomenologically 
uninteresting. Also, we would like to mention that in the 
present study, spin plays only a passive role, and we shall 
not discuss this particular aspect any further. 

Naumov [l2|, [ID, [3 [l5j] studied geometric phases for 
two and three flavor neutrino oscillations taking into ac- 
count the optic potentials [22j induced by coherent for- 
ward scattering of neutrinos against the background mat- 
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ter via SM interactions. The slowly changing parameters 
in the Hamiltonian were identified as a set of optic poten- 
tials q(t), which were connected to the refractive indices 
of neutrinos in a medium. For the naturally existing 
cyclic cases like spherically symmetric or sandwich-like 
density profiles, he found that the geometric (or topolog- 
ical) phase was zero for both two and three flavors due 
to only one of the optic potentials appearing in an essen- 
tial manner in the Hamiltonian. Note that the two terms 
"topological" and "geometric" were used interchangeably 
in Naumov's works. Here we will make a distinction be- 
tween the two terms. The topological phase refers to 
phase factors that are insensitive to small changes in 
the circuit, while geometric phases are sensitive to such 
changes. 

In a more recent paper, He et. al. [HI carried out a 
detailed study of the Berry phase in neutrino oscillations 
for both two and three flavors, active and sterile mixing, 
and with inclusion of nonstandard interactions. For the 
particular case of two flavor oscillations in matter, they 
claimed that the Berry phase can only appear if non- 
standard (R-parity violating supersymmetry) neutrino- 
matter interactions are taken into account. 

All the above papers [1 El El El EE El claim 
that the geometric phases do not appear in the oscil- 
lation probabilities for the case of two flavor neutrinos 
with CP conservation in vacuum or in matter as long as 
neutrino-matter interactions are standard i.e. coherent 
forward scattering is induced by charged current inter- 
action of electron neutrino (z/ e ) with electrons in matter. 
The above claims can be understood as the necessity of 
having at least two essential parameters in the Hamil- 
tonian to detect curvature. Because of the absence of 
flavor changing neutral currents in the SM, it turns out 
that for the case of ordinary electrically neutral matter, 
even though one has two varying parameters - electron 
number density (n e ) and neutron number density (n n ), 
only one of these will appear in an essential way in the 
Hamiltonian and hence the Berry's geometric phase is 
expected to be zero. The other parameter n n just adds 
a global phase to the time-evolved neutrino flavor state 
and hence does not affect oscillation. But also it is worth 
stressing that if both the conditions of having a nontriv- 
ial multidimensional parameter space as well as cyclic 
evolution of the states in parameter space were satisfied, 
the net geometric phase (resulting from the difference be- 
tween the geometric phases picked up by the individual 
mass eigenstates) would have appeared in the formulae 
for detection probability and hence been observable. 

Next, we will briefly review and summarize papers 
dealing with geometric phases that are generalizations 
of the Berry phase [20j, 123, |2J| in the context of neutri- 
nos [H EH- Such geometric effects can appear under 
less restrictive conditions than those required for Berry's 
version of the geometric phase. Infact such phases can 
appear even in situations where there are no parameters 
varying in the Hamiltonian and the evolution is not nec- 
essarily cyclic or unitary. Note, however, that in general 
the geometric phases appearing in transition amplitudes 



are global phases that do not have any observable conse- 
quences. To observe such a phase one needs a split-beam 
interference experiment in which a beam is spatially sep- 
arated into two parts that suffer different histories. Such 
an experiment is hard to design for neutrinos because 
they interact so weakly and are nearly impossible to de- 
flect or confine. This renders such phases uninteresting as 
they are not observable as far as neutrinos are concerned. 
Our aim here is to explore whether there are geometric 
effects that survive at the level of detection probabilities 
that are directly measurable quantities. 

Blasone et. al. [ItJ claimed that Berry's phase was 
present in the physics of neutrino oscillation in vacuum 
even for the two flavor CP conserving case. Their argu- 
ment is based on the fact that under Schrodinger evolu- 
tion, the pure flavor states come back to themselves after 
one period (T) of oscillation having acquired an overall 
phase. This overall phase was shown to be a sum of a 
pure dynamical phase and a part that depended on the 
mixing angle only and independent of energy and masses 
of the two mass states (hence, geometric). They called 
this extra phase the Berry phase. Note that this phase 
picked up by a neutrino flavor state arises purely due to 
Schrodinger evolution of the system giving a closed loop 
in the Hilbert space but not due to any slowly varying 
parameters leading to adiabatic evolution of the Hamil- 
tonian itself. Hence strictly speaking it is the Aharonov- 
Anandan cyclic phase [2(| that generalized Berry's adi- 
abatic phase to situations where the adiabaticity con- 
straint did not apply and only the cyclic condition is 
met. Also, we should note that since the phase obtained 
was a global phase at the amplitude level, it does not ap- 
pear in measurable quantities like neutrino appearance 
or survival probabilities as mentioned above. 

After Berry's [19| seminal paper on this subject, Ra- 
maseshan and Nityananda [25| pointed out that Berry's 
phase had a connection with the phase obtained by Pan- 
charatnam [23] in the fifties in his study on interference 
of polarized light. These insights were carried over to the 
ray space of quantum mechanics by Samuel and Bhan- 
dari |24j |. They showed that the two seemingly differ- 
ent geometric phases obtained by Berry and Pancharat- 
nam (appearing under different sets of conditions) could 
be described in a unified framework. They also pointed 
out that geometric phases are not restricted to unitary, 
cyclic and adiabatic evolution [l9| of a quantum system 
and can appear in an even more general context that al- 
lows for quantum collapses, which occur during measure- 
ments. Following this line of thought, Wang et. al. [H| 
extended the study of Blasone et. al. [17| to obtain non- 
cyclic geometric phases for two and three flavor neutrinos 
in vacuum. Their claim can be understood as follows. 
Consider the Schrodinger evolution of a quantum state 
over an arbitrary time period from t = to t. Now this 
open loop (noncyclic) Schrodinger evolution of a quan- 
tum state over a time r can be closed by a collapse of 
the time-evolved quantum state at r onto the original 
state at r = by the shorter geodesic curve joining the 
two states in the ray space j24l| . The phase associated 
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with the complex number (re zfi ) representing the inner 
product of the original state vector and the time-evolved 
state vector (with the dynamical phase removed) has a 
pure geometric origin. This noncyclic geometric phase 
was evaluated by Wang et. al. [l8| for both the two 
and three flavor cases. But, again note that this phase 
will be unobservable as it only appears at the level of 
amplitude. 

The main purpose of the present work is to estab- 
lish that Pancharatnam's phase does appear in detec- 
tion probabilities and hence is directly observable. For 
the simplest case of two flavors in vacuum or in con- 
stant density matter (restricting to SM interactions) with 
CP conservation, we obtain a Pancharatnam phase of 7r, 
and this leads to an elegant geometric interpretation of 
the neutrino oscillation formulae. We also make a direct 
connection of this phase with the Herzberg and Longuet- 
Higgins topological phase (26| in molecular physics. We 
show that the Pancharatnam phase of 7r remains even in 
the presence of slowly varying matter density and this 
can be ascribed to the topological nature of this phase. 
Inclusion of CP violation can change the topological na- 
ture of the phase and make it a path-dependent geometric 
phase. 

Although one should do a full three flavor analysis for 
a complete treatment, we work in an effective two fla- 
vor approximation that is fairly justified [I?], HH due to 
the smallness of 813, and hierarchy of mass splittings 
(l^mlx/^m^l << 1) and in addition on matter interac- 
tions being standard (50j . In many physical situations, 
observations depend on mainly one mixing and one mass 
squared splitting. Conventionally, 812 and 5m|i describe 
oscillations of solar neutrinos, while 023 and Srn^ 2 are 
used to describe atmospheric neutrinos. The mixing an- 
gle 0i3 gives small effects on both solar and atmospheric 
neutrinos. Working with only two flavors is of course 
advantageous as the results obtained are physically more 
transparent and can be visualized in analogous situations 
in optics and the Poincare sphere can be used as a cal- 
culational tool to study the system. 

For the ease of visualization of the phenomena of os- 
cillations, in the past several authors have discussed sim- 
ple pictorial depiction of neutrino oscillations in terms of 
precession of a (pseudo) spin vector in three-dimensional 
space in a variet y o f contexts for the case of two neu- 
trino flavors US, [M HI |M M, EH- Below we 
give a brief account of the papers dealing with geomet- 
ric representation of neutrino flavor oscillations. Harris 
and Stodolsky [29] addressed the question of a unified 
treatment of generic two-state systems (including parti- 
cle mixing involving two neutrino types) in media using 
density matrices. It was shown that the equation of mo- 
tion for the polarization vector represented the precession 
of polarization vector about a vector representing an ef- 
fective magnetic field (which could result from the mass 
terms in vacuum or matter terms) . Kim et. al. [30( dis- 
cussed the analogy of solar-neutrino oscillations with that 
of precession of electron spin in a time-dependent mag- 
netic field. They applied this picture in the limit of adia- 



batic approximation. Stodolsky [3l| described the evolu- 
tion of a statistical ensemble (neutrinos from supernovae 
or in the early Universe) applying the density matrix ap- 
proach [2t| and showed that oscillations in presence of 
mixing and matter interactions in a thermal environment 
could be viewed in terms of precession. Kim et. al. [32j 
derived the geometric picture for two and three flavor 
neutrinos and applied it to nonadiabatic as well as adia- 
batic cases. Thomson and McKellar [33| treated the case 
of neutrino background giving rise to nonlinear feedback 
terms in the equation of motion for polarization vectors 
and gave a pictorial representation for the same. En- 
qvist et. al. [34[ describe visualization of oscillations of 
a thermal neutrino ensemble of the early Universe. The 
geometrical representation in wave packet treatment of 
oscillations was discussed by Giunti et. al. [35}. As in 
optics, the Poincare sphere is a convenient tool for visu- 
alizations and calculations pertaining to neutrino oscilla- 
tions, particularly in looking for geometric effects. 

This paper is organized as follows. In Sec. [HI we- de- 
velop an analogy between the neutrino flavor states and 
polarized states in optics since such a mapping allows 
for a convenient visualization of geometric effects. We 
then go on to show in Sec. IIIII that the Pancharatnam 
phase does appear in the detection probabilities of neu- 
trino species in the two flavor neutrino system in vacuum 
and also in matter. We conclude with a discussion of our 
key result and future directions in Sec. HVl Throughout 
we set h = c = I. 



II. CORRESPONDENCE BETWEEN TWO 
FLAVOR NEUTRINOS AND POLARIZATION 
STATES IN OPTICS 

Since the concept of Pancharatnam's phase was devel- 
oped in the context of optics, it is worthwhile to first 
develop a correspondence between the mathematics of 
two flavor neutrino states and polarization states in op- 
tics. Let us first recall the conditions under which the 
two flavor neutrinos and polarization states in optics can 
be analyzed within an unified framework. 

A. Two flavor neutrinos 

In the ultra-relativistic limit, the Dirac equation for 
two flavor neutrinos (antineutrinos) can be reduced to 
a Schrodinger form [22 . [37} written in terms of a two- 
component vector of positive (negative) energy probabil- 
ity amplitude. This is analogous to Maxwell's equations 
reducing to the linear Schrodinger form for the polariza- 
tion states in optics in the paraxial limit [38] . 

The two neutrino flavor states can be mapped to a 
two-level quantum system with distinct energy eigen- 
values, Ei ~ p + mf /2p in the ultrarelativistic limit 
along with the assumption of equal fixed momenta (or 
energy) [H, In the presence of matter, the relativis- 
tic dispersion relation Ei — f(j>, mi) gets modified due to 
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the neutrino-matter interactions (in an electrically neu- 
tral homogeneous medium) leading to 

/ mf + mf. V c Jr \ 

T^\f(ujsm20 f + (V c - wcos26) 2 , (1) 

where u> = 8m 2 /2p with mass splitting 8m 2 — m 2 — m 2 
and p ~ -E being the fixed momentum (energy) of the 
neutrino. 6 is the mixing angle in vacuum. Vc — 
V2G F n e = 7.6 x lCT 14 y e p eV and V N = -V2G F n n /2 = 
—3.8 x 10~ 14 Y n p eV are the respective effective poten- 
tials due to coherent forward scattering of neutrinos with 
electrons (via charged current interactions) and neutrons 
(via neutral current interactions). Gp — 1.16637 x 
1CP 5 GeV -2 parameterizes the weak interaction strength 
(Fermi constant). Vc and Vn depend on the electron 
(n e ) and neutron (n n ) number densities (in units of 
cm~ 3 ). n e / n — pY e / n NAvo, where p is the mass den- 
sity in g cm -3 , Y e / n is the relative electron (neutron) 
number density and its value is roughly ~ 0.5 for Earth 
matter, and Navo is the Avogadro's number. Setting 
Vc — Vn — 0, we recover the vacuum case. 

Note the fact that although there are two densities n e 
and n n appearing in the eigenvalues, it is only n e that 
appears in a nontrivial way (through Vc) in the flavor 
Hamiltonian, 

/ m? + m 2 . Vn \ 

1 / V c - u cos 26 u sin 26 \ , , 
2^ wsin26 -(Vc - w cos26) J ' > 

The above Hamiltonian (Eq. [2]) also describes an inhomo- 
geneous medium provided the scale of variation of mat- 
ter induced potential Vc is slow compared to the scale 
of the order of h/(E + — EJf), hence ensuring no tran- 
sitions between the mass eigenstates. This defines the 
adiabaticity condition (36l . [33 |. As neutrinos traverse a 
density gradient, at a particular value of n e the diagonal 
elements of H„ can vanish causing an interchange of fla- 
vors irrespective of the value of the vacuum mixing angle 
6. This phenomenon of resonant conversion in matter 
is known as the Mikheyev-Smirnov-Wolfenstein (MSW) 
effect 0, El|. 

The off-diagonal form of the Hamiltonian in flavor ba- 
sis (both in vacuum and matter) leads to flavor oscil- 
lations of neutrinos, which is the only mechanism that 
mixes the neutrinos of different generations or flavors 
while preserving the lepton number (note that the ab- 
sence of flavor changing neutral currents prevents any 
flavor change within the SM). Also note that the matter 
term appears in diagonal elements only so in the absence 
of vacuum mixing, neutrinos of different flavors cannot 
mix. The term proportional to the identity gives an over- 
all phase to each of the mass eigenstates and hence does 
not affect oscillations. This corresponds to the gauge 
freedom of any state of a two- level quantum system [3j . 



In the next subsection, we describe the polarized states 
in optics in the language of quantum mechanics. 



B. Polarized states in optics 

Polarization optics is mathematically identical to the 
evolution of a two state quantum system. In a helicity 
basis for polarized light, we can write | 9t) and | £) rep- 
resenting right and left circular polarizations. A general 
polarized light beam | ^ ) can then be expanded in this 
basis as | * ) = a\ 9t> +/3| £) where \a\ 2 + \(3\ 2 = N, the 
intensity of the beam of polarized light. We can param- 
eterize an arbitrary state of polarized light by 

\W) - ViVexp^l ^ sm (0/2)exp(#/2) J' [6) 

where N is the total intensity, which is normalized to 
unity, and the angles 9 and <j) (where < 9 < -k and 
< (f> < 2ir) describe the state of polarization of the 
beam, represented on the two-dimensional unit sphere 
(S 2 ) called the Poincare sphere. Orthogonal polarization 
states are antipodal points of the sphere, ij is the overall 
phase of the beam. The states on the sphere are de- 
fined modulo this overall phase of r\ and represent the 
ray space 42]. The north pole (9 = 0) represents right 
circular light and the south pole (9 — n) represents left 
circular light. States on the equator (0 — 7r/2) repre- 
sent linear polarizations. Any other point on the sur- 
face of the sphere represents elliptic polarization. The 
Poincare sphere is a useful device to visualize the changes 
in the state of polarization of a light beam traversing 
through a medium. 

The mapping between the polarized states and a two- 
level quantum system originates from the following fact. 
Neglecting absorption effects [51(, the effect of different 
media can be encoded in terms of 2 x 2 Hermitian ma- 
trix (Hamiltonian). The time evolution of optical states 
in a medium is governed by a Schrodinger-like equation 
with the medium represented by the most general form 
of Hamiltonian for a two-level system given by 

H = Aa x + Ba v + Ca z + VI , (4) 

where, the coefficients of the three traceless Pauli ma- 
trices, A, B and C are responsible for generating rota- 
tions of incident optical states about x,y, z axes on the 
Poincare sphere. T> just adds an overall phase that can 
be absorbed in a redefinition of the state. Hence given 
an arbitrary medium, it can be represented by a Hamil- 
tonian as mentioned above, and the eigenstates of the 
Hamiltonian represent those optical states that do not 
suffer any change (when incident on such a medium) in 
their state of polarization except for picking up an overall 
phase shift. The polarization of any other state (other 
than the eigenstates) incident on this medium will un- 
dergo a periodic change. On the Poincare sphere this can 
be visualized as a rotation of the incident state vector 
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FIG. 1: Neutrino states on the Poincare sphere. The flavor 
states | v a ) and | vp ) are the two antipodal points on the z 
axis while | #, ± ) correspond to the mass (energy) eigenstates 
lying on an axis making an angle i9 with respect to the z axis. 



about the axis denned by a line joining the two eigen- 
states of the Hamiltonian. Mathematically, these uni- 
tary rotations on the Poincare sphere are generated by 
e~ %m . This is identical to unitary time evolution gen- 
erated by the Hamiltonian of the quantum states in the 
Hilbert space. The quantum-mechanical analogue of the 
Poincare sphere is the Bloch sphere, which geometrically 
represents the space of pure states of a two-level quantum 
system. 

Nonvanishing values of A, B,C simultaneously param- 
eterize the effect of an elliptically birefringent medium. 
Circular (linear) birefringence are special cases where 
the conditions A,B = and C,X> ^ (B,C = and 
A, T> 7^ 0) are satisfied. 



rotating the north pole into the south pole. 

Ignoring the term proportional to the Identity, the neu- 
trino Hamiltonian (Eq. [2|) both in vacuum or matter can 
be recast in exactly the same form given by (see Eq. U) 



H <> = 2 [(sini?)^ - (cosi?)^ 



(5) 



where ui = Sm 2 /2p and the mixing angle is replaced 
by -0/2 [H^]. Comparing the two Hamiltonians (Eq. 0] 
and Eq. [5]) we see that the neutrino Hamiltonian repre- 
sents a medium with elliptic birefringence. And neutrino 
oscillations can be viewed as the neutrino flavor state 
precessing (3(| about the line joining the mass eigenstates 
(analogous to elliptic axis) induced by the time-evolution 
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on the Poincare sphere. In the language 



of neutrino optics, both vacuum and matter exhibit el- 
liptic birefringence property with different elliptic axes. 

The absence of flavor changing neutral currents in the 
SM gives rise to a real form of the Hamiltonian [B = 0), 
and it corresponds to a CP-conserving situation. The 
eigenvectors (also called mass eigenstates) of Eq. [5] are 
given by 



/cosW2)\ and |^ x f-sinW2)\ 
" ^sin(tf/2) ) 1 ' >~\ cos(tf/2) ) ■ 

(6) 

Note that states | + ) and | — ) are orthogonal an- 
tipodal points on the Poincare sphere which always lie 
on the great circle formed by the intersection of the x — z 
plane with the Poincare sphere. Mass eigenstates lying 
outside the x—z plane imply CP violation. This fact has 
very interesting consequences for the physics of geometric 
phases in CP nonconserving situations [43| . 



III. PANCHARATNAM'S PHASE IN THE TWO 
FLAVOR NEUTRINO SYSTEM 



C. Neutrinos and optics analogy 

We can now describe the isomorphism between neu- 
trino states and polarized states in optics. The complete 
set of states for two flavor neutrino system can be repre- 
sented on the Poincare sphere just like the optical states 
as depicted in Fig [TJ For convenience we define a new 
coordinate d, which goes from — » 2-7T as we traverse 
the unit great circle in the x — z plane. In terms of the 
old coordinates, the points 0, <fi — are now labeled by 
■0 = and the points 0, <f> = ir are labeled by ^ = 2tt — 0. 
If we assume that the flavor states are the north and 
south poles of the Poincare sphere, then the mass eigen- 
states are represented by the two antipodal points ly- 
ing on an axis making an angle 20 = i0 with respect to 
the polar axis. States on the equator coincide with the 
mass eigenstates for the special case of maximal mixing 
(0 = -0/2 = 7r/4) which corresponds to complete flavor 
conversion (MSW effect). Geometrically, the MSW ef- 
fect can be viewed as rotation about an equatorial axis, 



The Pancharatnam phase :- We give a brief introduc- 
tion to the idea of Pancharatnam's phase in quantum- 
mechanical language along the lines of Ref . [24|, |42j, |4J] . 
Given any two nonorthogonal states | 21 ) and | 03 ) in the 
Hilbert space describing a system, a notion of geometric 
parallelism between the two states can be drawn from 
the inner product ( 21 1 58 ) . The two states are said to be 
parallel (in phase) if ( 21 1 58 ) is real and positive, which 
defines the Pancharatnam connection ( or rule ). Geomet- 
rically, it implies that the norm of the vector sum of 
the two states ||(|2l) + |58))|| 2 = (2t|2l) + (5B|<B) + 
2| ( 21 1 58 ) | cos(ph( 21 1 58 )) is maximum. Physically, it im- 
plies that if we let the two states interfere with each other 
the resulting state will have maximum probability (inten- 
sity) . Note that if | 21 ) is in phase with 1 58 ) , and 1 58 ) is 
in phase with | £ ) , then | £ ) is not necessarily in phase 
with state | 21 ) . The phase difference between the states 
| £ ) and | 21 ) is the Pancharatnam phase, and it is equal 
to half the solid angle Vl subtended by the geodesic trian- 
gle 21, 58, € on the Poincare sphere for a two- level system 
at its center. In general, for an n-level system, the space 
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of states is given by CP n_1 (CP stands for complex pro- 
jective) which reduces to the Poincare sphere (S 2 ) for a 
two- level system (n = 2). Nonintegrability of Pancharat- 
nam's connection follows from the nontransitivity of the 
rule. 

Pancharatnam's phase reflects the curvature of pro- 
jective Hilbert space (ray space) and is independent of 
any parameterization or slow variation. Thus it can 
also appear in situations where the Hamiltonian is con- 
stant in time. All one needs is that the state has 
a nontrivial trajectory on the Poincare sphere. This 
condition is met naturally for neutrinos since they 
are produced and detected as flavor states (which are 
not the stationary mass eigenstates) and hence they 
automatically explore the curvature of the ray space 
(Poincare sphere) under the Schrodinger time evolution. 
Furthermore, note the fact that Schrodinger evolution 
(possibly) interrupted by measurements can lead to Pan- 
charatnam's phase. If we take any state and subject 
it to multiple quantum collapses (such that consecutive 
collapses are between nonorthogonal states) and bring 
it back to itself, then the resulting state is given by 
|2l)(2l|<£)(£|Q3)(58|2l), where the phase of the com- 
plex number (2l|£)(e:|<8)(<8|2l) is given by ft/2. 

The Herzberg and Longuet-Higgins phase and CP- 
conserving neutrino Hamiltonian :- Let us reexamine 
the form of the neutrino Hamiltonian given by Eq. [5] and 
the eigenvectors given by Eq. [S] Note that the eigenvec- 
tors depend only on a single parameter § and satisfy 

|tf,±) = T|tf + 7r,T) = -|0 + 2tt,±) 

= ±|tf + 3tt,t} = |i? + 4tt,±) . (7) 

The minus sign picked up by both the mass eigenstates 
as we change •& from — > 2tt is precisely the Herzberg 
and Longuet-Higgins phase [26l . |45| | of tt, which was first 
obtained in the context of molecular physics in 1963. So, 
we note that just by looking at the form of the Hamilto- 
nian for neutrino system, we should expect the Herzberg 
and Longuet-Higgins phase to appear. Also, note that 
the space of rays for the real neutrino Hamiltonian is 
the great circle (S 1 ) lying on the x — z plane of the 
Poincare sphere (Fig.[T]) and global structure of the eigen- 
vectors is a Mobius band. The variation of "d results in 
parallel transport of the mass eigenstates (with dynam- 
ical phase removed) following the parallel transport rule 
along $ , 

SM*fI^I*f> = 0. (8) 

This parallel transport rule (formally referred to as nat- 
ural connection) has an anholonomy defined on the 
Mobius band and this leads to the topological phase of 
tt. The topological phase factor 6 depends on the vector 
potential A§ given by 

Q = <f Afid-d = (f Qm( $ T \ -^-\d T ) d-d . (9) 



This vector potential A$ is nonintegrable, and this is the 
anholonomy of the connection. Physically, this corre- 
sponds to half a unit of magnetic flux piercing the origin 
of the x — z plane, encircling which leads to this topolog- 
ical phase. And, the origin of the circle is connected to 
the null Hamiltonian (i.e. all elements are zero), which 
corresponds to the degeneracy point. 

Naively speaking, one would think that this phase will 
be impossible to access for neutrinos because we do not 
have a handle on the mixing angle i?/2 to be varied in 
a controlled way from = — * 2ir. The key point to 
understand here is the fact that as long as we carry out 
a quantum evolution of a state in a closed loop enclos- 
ing the point of singularity (degeneracy point, origin of 
the Poincare sphere), which can be achieved either via 
adiabatic variation of i? or via Schrodinger evolution in- 
terrupted by collapses, one will always get this phase. 
However, note that in the former case, the amplitude of 
the initial state undergoing evolution does not change 
but in the latter case, it diminishes. In what follows, 
we will show that the transition probability for neutri- 
nos actually does carry imprints of such a topological 
phase, which can be explicitly derived using Pancharat- 
nam's prescription. We then show that the phase of tt 
actually appears there and is in fact observed by all the 
experiments carried out so far. 

The topological phase in two flavor neutrino oscilla- 
tions (invoking collapses and adiabatic evolution) :- In 
what follows, we consider the most general situation, i.e. 
neutrinos are traversing through matter with slowly vary- 
ing density (i.e. & is a slowly varying parameter changing 
from $i to f^)- Vacuum or constant density matter will 
be special cases where d is a constant. 

In order to see the effect of geometric phases, usually 
one performs a split-beam experiment. In the case of op- 
tics, one separates a beam into two parts in space and 
each part traverses a different path. Finally the beams 
are recombined to observe the relative phase shift as they 
interfere. In optics, the reflective and refractive prop- 
erty of the medium is exploited to make devices like mir- 
rors and lenses, which facilitates designing of such exper- 
iments in the laboratory. In the case of neutrinos, such 
a procedure is not possible owing to the fact that the re- 
fractive index is extremely small (n re f r — 1 ~ 10~ 19 for 
neutrinos of energy 1 MeV in ordinary matter). Treat- 
ing the Sun (with density p = 150 g cm -3 in the core) 
as a spherical lens for a neutrino beam of energy 10 MeV 
passing through it, one gets the focal length to be around 
lO 18 i?0 [1^, which is about 10 5 times the size of our 
galaxy. Spatially split-beam interference experiments 
with neutrinos are clearly impossible. However, the fact 
that neutrinos are produced and detected as flavor states 
allows us to think of the time evolution of neutrinos as 
a split-beam experiment in energy space as depicted in 
Fig.[H 

Let us consider a neutrino created as a flavor state 
| u a ) (for example, neutrinos produced inside the Sun 
are predominantly in the electron neutrino flavor state, 
| v e )) and detected as another flavor state, \vp) (\vp) 
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FIG. 2: Schematic of a split-beam experiment for neutrinos in 
energy space. | v a ) and 1 vp } are the two flavor states, while 
| ± ) and also | $2, ± ) correspond to two sets of mass (en- 
ergy) eigenstates. [ <&\, ±) are adiabatically evolved to states 
| #2, i), respectively (upon removing the dynamical phase). 



can either be a | ^ e ), i.e. survival of the same electron 
neutrino flavor or a | ), i.e. appearance of muon neu- 
trino flavor), then 

I f a ) =V a+ \<&!, + ) + , (10) 

where |i?i,±) are the eigenstates of H„(#i). Now we 
consider an adiabatic evolution of the mass eigenstates 
from | #i, ± ) to | $2, ± ) due to a slow enough variation of 
background density such that no mixing between the two 
eigenstates is ensured under time evolution, and | ± ) 
evolves to 

|i?i,±) e- lp ±|^ 2 ,±) with 

V± = ±i /" ^(w sin i?) 2 + ( Vc - w cos i?) 2 (ft' 

+ {>(l + ^i + v f + VN y An) 

as the dynamical phases, relevant both for the vacuum 
case (Vc = Vn = 0) and in the presence of varying mat- 
ter density profile and t is the time of flight of the neu- 
trino. The quantities that depend on time (or distance) 
are Vc and Vn defined earlier (see Eq. [I}. Note that 
the states | ± ) are | "d 2l ± ) are connected via parallel 
transport rule (Eq. [5]) on the Poincare sphere. The two 
time-evolved states e~ lV± | $2, ±) are finally recombined 
to form a flavor state at the detector. 

In order to see this explicitly, let us proceed as follows: 
The amplitude for the transition between states v a — > vp 
is given by 

A(v a -> up) = (vp \U\v a ) , (12) 

where IA is the unitary evolution operator given by 

U = e- w +\$ 2l + ){d ll + \+e- w -\d 2l -){d ll -\ . 

(13) 



Inserting two complete sets of states in the amplitude, 

+ 

A{v a ^v f} )= (vp\#2,i)(#2,i\U\AJ)(#i,j\v a ) 

= (vp\$2, + ){'&2, + \U\& 1 , + ){§ 1 , + \v a ) 

+ („ \'& 2 ,-){'& 2 ,-\U\d 1 ,-){§ 1 ,-\v a ) . (14) 

Note that the cross terms do not contribute in the adia- 
batic limit. Upon substituting Eq. [T3]in Eq. [TJJ we get 

A{v a ^vp) = e- i75 +(^|tf 2 , + )(tfi, + K) 

+ e- w -(vp\d 2 ,-)(d x ,-\v a ) . (15) 

Then the probability for flavor transition v a — > vp is 
given by 

V{v a —> vp) = \A{v a -> vp)\ 2 

= (v a \0u + )(Qn, + \vi})(vfi\$2, + )(di, + \v a ) 
+ {v a \d ll -){d 2l -\ V p){vp\d 2 ,-){d l ,-\v a ) 

+ [(v a \4i,- )e lV - (02,-\i'i3){vfi\02, + )e- W+ 
(4i, + \v a )+c.c.]. (16) 

The cross term term in Eq. [TO] is related to the inter- 
ference term resulting from the two path interferometer 
depicted in Fig. O Upon dropping the dynamical phase, 
we have ( v a | - )(i? 2 , - | vp )( vp \ <d 2 , + )( + | v a ) 
which can be viewed as a series of closed loop quantum 
collapses with intermediate adiabatic evolutions given 
by \v a ) - + ) -> |0 2 , + ) - |^) -» |tf 2 ,-) -> 
I — ) — > I f a ) that essentially covers a great circle in 
the x — z plane as is shown in Fig. [21a). This closed 
trajectory subtends a solid angle of SI = 2-tt at the center 
of the great circle. Hence without any further calcula- 
tion, we can immediately predict that the phase of the 
interference term will be ir (half the solid angle) due to 
Pancharatnam's prescription. On the circle, each of the 
individual collapse processes which essentially projects a 
state with given angle i9 to another state with different 
angle can be thought of as an infinite series of infinites- 
imally close collapses between states defined as 1 1) ) and 
I fl + 6$ ) as far as geometric phases are concerned. The 
entire closed loop of collapses with intermediate adiabatic 
evolutions mentioned above can be viewed as a smooth 
variation of i9 from — > 2tt in the limit 8§ — > hence 
making a direct connection to the Herzberg and Longuet- 
Higgins phase mentioned above. Nonetheless, we must 
note that the evolution of a state is unitary under in- 
finitesimal collapses (Stf — > limit) while it is nonunitary 
under finite collapses leading to a loss in intensity (prob- 
ability). But the geometric phase of the evolving state 
remains unaltered for the two cases mentioned above. 

For the case when a — (3, i.e. survival probability, 
it is easy to see that the collapses do not form a closed 
loop enclosing the origin and therefore the interference 
term will not pick up any phase. This case is depicted in 
Fig.^b). 

In a simpler situation when ■d does not change, i.e. the 
case of vacuum or constant density matter, the number 
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| Vac ) \v a ) 




FIG. 3: Two representative cases depicting the collapse processes (dashed red lines) with intermediate adiabatic evo- 
lutions upon removing the dynamical phase (dotted blue lines) on the great circle (S 1 ) arising due to the cross term 
( Va. 1 #1, — ) ( $2, — | Vp ) { vp | $2, + ) ( 0i, + | v a ) in the probability. The initial flavor state | v a ) is on the positive z axis, while 
the final flavor state \ vp) is not necessarily its antipodal point. The two sets of mass eigenstates are antipodal points on 
two axes making angles #i and #2 respectively with respect to the 2 axis. Case (a) corresponds to appearance probability 
\P(y a — * Vp)] for which we get a cyclic loop in $ space, (b) The collapse processes for survival probability \P(y a — > v a )] do not 
enclose any loop. 



of states will be fewer (in the absence of variation of den- 
sity, I $1, ± ) is the same as | #2, ± )) and the collapses are 
given by | v a ) -> | <? 1} + ) -> | vp ) -> | - ) -> | v a ). As 
long as the collapses lead to closed loop encircling the ori- 
gin, we will obtain this topological phase. So this phase 
of 7r appears whether we consider vacuum and/or ordi- 
nary matter with constant density or with slowly chang- 
ing (but noncyclic) electron number density. This is due 
to the topological character of this phase, which will be 
preserved as long as we have CP-conserving (real) Hamil- 
tonian and states are always lying on a great circle in the 
x — z plane in the Poincare sphere. 

Next we write down an explicit expression for the ob- 
servable quantities, i.e. appearance and survival prob- 
abilities for two neutrino flavors. Using the general ex- 
pression obtained in Eq. 1161 the appearance probability 
for transition v e — ► is given by [531 ] 

v{ Ve ^^) = u: + (e 1 )u A1+ (e 2 )u* + (e 2 )u e+ (e 1 ) 
+ u:_(e 1 )u Al _(e 2 )u*_(e 2 )u e _(e 1 ) 
+ [u:_(e 1 ) e lI, -u Al _(e 2 )u; i+ (e 2 ) e - lI, +u e+ (e 1 ) + c.c.]. 

(17) 

Note that the matrix U(0) is the lepton mixing matrix 
(defined in a basis where the charged lepton mass matrix 
is diagonal) . It is also referred to as the Pontecorvo-Maki- 
Nakagawa-Sakata (PMNS) matrix (4|| [4lJ and connects 
the flavor states to the mass eigenstates. For the 2x2 
case, it is a real orthogonal rotation matrix given by 

U(6)= f C ° s6 n Sin ^) . (18) 
v ' \— sin0 cos By v 7 

Substituting the elements of U(0) we get 

Viv e ^v il ) = cos 2 9i sin 2 2 + sin 2 0i cos 2 9 2 
+ [2cos(£>+ -£>_)](- sin0i)cos© 2 sin 2 cos0i (19) 



We note that there are four inner products appearing in 
the interference term in the final expression for the prob- 
ability out of which the first three inner products, viz., 

{0l, + \v e ) - U e+ (0!) = COS©! > 0, (l/„|02, + ) = 

U*+(0 2 ) = sin 2 > and {^ 2 ,-|^) = U M _(0 2 ) - 
cos 02 > clearly implying that these states are mu- 
tually parallel to each other in pairs according to Pan- 
charatnam's rule, which is to have the inner product 
of any two states real and positive, while the last one, 
{v e \di,-) = U*_(0i) = -sin©i < by Pancharat- 
nam's rule has \ v e ) antiparallel to — ), since the 
physically allowed values for the mixing angles ©i and 
02 are within the interval [0,7r/2] for Sm 2 > [3^| (On 
the Poincare sphere, the corresponding $1 and $ 2 can 
take values between [0, w]). The minus sign appearing in 
the interference term is thus the Pancharatnam's phase 
of 7r appearing in the neutrino oscillation formula (see 
Fig. [Ha)). 

If in a hypothetical situation, for some range of param- 
eters 0i and 02, the first three of the inner products are 
real and negative (i.e. the states are aligned antiparal- 
lel to each other or completely out of phase), while the 
fourth inner product is real and positive (the states are 
in phase) then also we will have this minus sign. The 
nontransitivity also holds here leading to the non-trivial 
topological phase of %. This situation where the inner 
product becomes real and negative defines an "antipar- 
allel" rule (in the same spirit in which Pancharatnam 
defined his rule of two states being "in phase or paral- 
lel" ) would correspond to the norm of the vector sum of 
the two states being at its minimum value. Physically, 
this implies the interference of the two given states would 
be destructive and the resulting state will have minimum 
intensity or a dark fringe in optics. 

The existence of Pancharatnam's phase of it can be 
simply connected to the fact that the mixing matrix U(0) 
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matrix for two flavors is an orthogonal rotation matrix 
parameterized by the mixing angle O of which one ele- 
ment has a negative sign. Thus, this phase is built into 
the structure of U(O) matrix. 

The survival probability is given by 

V{v e — ► v e ) = cos 2 0i cos 2 2 + sin 2 0i sin 2 6 2 
+ [2cos(X>+ -£>_)] sin 0i cos 2 sin 2 cos 0i . (20) 

Note that in the case of survival probability, the cross 
term does not pick up any nonzero topological phase, and 
geometrically this is exactly what we had expected from 
FigCJb). The loop in z9-space is open in this case, and 
this is what leads to this result. The topological phase of 
the interference term in survival probability is zero, while 
it is 7r in the case of the appearance probability, and this 
fact is in accord with unitarity. 

The above expressions (E qs. IT9l and |20|) reduce to the 
standard results [27|, |36|, |39j, |48( for vacuum if we substi- 
tute 0i = 2 = 0, 

_, . . 9 . o Sm 2 l 

r{y G — > Vu) — sin 20 sin and 

AE 

Sin 2 1 

V{u e ~^v e ) = 1 - sin 2 20 sin 2 , (21) 

where in the ultrarelativistic limit, we can use t ~ I and 
p ~ E leading to T>± = ±Sm 2 l/2E (see Eq. [TT]) for the 
vacuum case (Vc — Vn = 0). In constant density mat- 
ter, the quantities and Sm 2 in Eq. [5T]are replaced by 
their respective renormalized values in matter, m and 
(Sm 2 ) m but the form of the expression will remain the 
same. Hence our result is consistent with the standard 
neutrino oscillation formulation, and it provides a clear 
geometric interpretation of the phenomenon of neutrino 
oscillations. More precisely, the standard result for neu- 
trino oscillations is in fact a realization of the Pancharat- 
nam topological phase. 

IV. DISCUSSION 

As mentioned in the introduction, the existing work 
on the subject of geometric phases in neutrino oscilla- 
tions led to the widespread belief that the two flavor 
neutrino oscillation formulae in CP conserving situations 
were devoid of any geometric or topological phase com- 
ponent. Appearance of the cyclic Berry phase was dis- 
missed on the grounds of not having any time-varying 
parameter in vacuum and having only one essential pa- 
rameter (thereby enclosing no area) in the case of normal 
matter d, [12, 0, [3 [H, [l6| . Concerning the appearance 
of the general geometric phase in the two flavor neutrino 
case for propagation in vacuum, there are claims report- 
ing its appearance [13, [HI]. But, it should be noted that 
such terms appeared only at amplitude level and as ar- 
gued earlier, a phase appearing in the amplitude can be 
observed only via a split-beam experiment, which is not 
feasible to design in the case of neutrinos. 



In this paper, we have examined the minimal case 
of two flavor neutrino oscillations and CP conservation. 
Contrary to all existing claims in the literature concern- 
ing the geometric or topological phase in two flavor neu- 
trino oscillation probabilities, our study provides the first 
clear prediction that a topological phase of n exists at the 
probability level even in the minimal case of CP con- 
servation. We show that it is inherently present in the 
physics of neutrino oscillations via the structure of the 
PMNS neutrino mixing matrix. This existence of this 
topological phase is linked to the presence of a flux line 
of strength ir at the origin of ray space, which is con- 
nected to the degeneracy point associated with the null 
Hamiltonian. 

Pancharatnam's idea is quite useful in terms of pre- 
dictive power as it allows for a clear visualization of the 
appearance of such a phase due to geometric effects with- 
out doing any algebra. Our prescription is general as it 
contains effects due to collapses and also due to adiabatic 
evolution. In the absence of either of these, one would 
get the same phase. So no matter what the details are, 
as long as the singular (degeneracy) point is enclosed by 
a cyclic loop (in the space of rays) as i9 is varied from 
— > 2ir, we will get this phase, and this is due to its 
topological robustness. The adiabatic and collapse pro- 
cesses both conspire in such a fashion that the net phase 
would always be n. This does not happen for geometric 
phases. 

The topological phase obtained in this paper is a con- 
sequence of anholonomy, which can arise in situations 
even when there is no curvature. The most striking ex- 
ample of this is the Aharonov-Bohm effect [4|| . To expe- 
rience the effect of anholonomy, the main requirement is 
to encircle the singular point, this fact was exploited by 
Herzberg and Longuet-Higgins in pointing out the topo- 
logical phase in molecular physics. On the other hand, for 
Berry's phase to appear, a net curvature is a must which 
is fulfilled by having at least two essential parameters in 
the Hamiltonian varying cyclically. This is an important 
distinction between the geometric phases as obtained by 
Herzberg and Longuet-Higgins and by Berry. 

If we consider mixed flavor states instead of the 
pure flavor states, there will be a greater number of phys- 
ical situations (or, possible diagrams for the interference 
term like the ones shown in Fig. [3] for pure flavor states) 
that can be explored to see if one encircles the singular 
point or not. A mixed state corresponds to a general 
point on the surface of the Poincare sphere like an el- 
liptically polarized state in optics. If the mixed states 
are such that they lie on the x — z plane, it will al- 
ways lead to the same quantized topological phase of 7r. 
But, for a general mixed state lying anywhere else on the 
Poincare sphere, the phase will be geometric in nature. 

It might be a nontrivial task to extend our geometri- 
cal interpretation to the case of three neutrinos flavors 
because it will involve a higher dimensional sphere (the 
ray space is CV 2 for the three level quantum system). 

It is natural to ask what happens when we invoke CP 
violation. In vacuum, CP violation cannot be induced 
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in the two flavor case as a consequence of CPT invari- 
ance and unitarity [28| . However, matter with constant 
or varying density can induce CP violation via the co- 
herent forward scattering of neutrinos with background 
matter. If we introduce CP violation induced by back- 
ground matter with constant density [28| , we still expect 
to get the same phase of tt as we have two pairs of or- 
thogonal states that will always lie on a great circle. If 
the density is varying slowly (adiabatic condition holds), 
then the intermediate states (connected by adiabatic evo- 
lution) will be lifted from the great circle, hence resulting 
in a path-dependent solid angle, and the phase will be ge- 



ometric [43l | . 
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